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We investigate the role of nonperturbative, bubble-like inhomogeneities on the decay rate of false- 
vacuum states in two and three-dimensional scalar field theories. The inhomogeneities are induced 
by setting up large-amplitude oscillations of the field about the false vacuum as, for example, after 
a rapid quench or in certain models of cosmological infiation. We show that, for a wide range of 
parameters, the presence of large-amplitude bubble-like inhomogeneities greatly accelerates the de- 
cay rate, changing it from the well-known exponential suppression of homogeneous nucleation to a 
power- law suppression. It is argued that this fast, power-law vacuum decay - known as resonant nu- 
cleation - is promoted by the presence of long-lived oscillons among the nonperturbative fluctuations 
about the false vacuum. A phase diagram is obtained distinguishing three possible mechanisms for 
vacuum decay: homogeneous nucleation, resonant nucleation, and cross-over. Possible applications 
are briefly discussed. 



I. INTRODUCTION 

Since the seminal results by Coleman which can 
be viewed as an extension of Langer's theory of Homo- 
geneous Nucleation (HN) Q in condensed matter sys- 
tems to relativistic field theories, there has been a large 
amount of work dedicated to vacuum decay at both zero 
and, inspired by Linde's work finite temperatures [1]. 
Several textbooks and review articles describe the topic 
in detail [1]. 

In high energy physics, the interest in false vacuum 
decay comes from the fact that models describing the 
fundamental interactions of matter fields often possess 
metastable states. For instance, it is still unclear if the 
quark-hadron transition is or not at least weakly first or- 
der @. Within models of electroweak symmetry break- 
ing, several extensions of the Standard Model, supersym- 
metric or not, support metastable states Q- It is hoped 
by many that results from the Large Hadron Collider will 
shed light on this issue, as they may reveal the fundamen- 
tal mechanism of mass generation. As we move toward 
the early universe, many models of inflation, including 
the original old inflation and many others, make use of 
potentials with metastable states The same is true 
of grand unified theories. It is thus of great interest to 
examine under what conditions the predictions from HN 
theory, which are widely used in the literature, can be 
trusted. In this work, we will explore the mechanisms 
by which a false vacuum can decay. As we hope to con- 
vince the reader, the effective decay rate is sensitive to 
the properties of the initial state: the mechanism of false 
vacuum decay reflects its previous history. Thus, hav- 
ing information about the decay mechanism may help us 



reconstruct the conditions prevalent at an earlier epoch, 
when the system under study was still in its metastable 
state. 

This paper is organized as follows: in the next section, 
we introduce our model and review basic results from 
false vacuum decay theory. We also briefly review the 
properties of oscillons, long-lived, time-dependent non- 
perturbative field configurations, as they play a key role 
in the present work. We complete the section describ- 
ing the Hartree effective potential and the details of our 
3d lattice implementation. In section HI, we explore the 
different mechanisms of vacuum decay, emphasizing the 
departures from HN. We describe in detail how a power- 
law decay rate is operative for a wide range of parameters 
controlling the properties of the initial state and the field 
interactions. We construct a phase diagram encompass- 
ing the three roads toward vacuum decay: HN, resonant 
nucleation, characterized by power-law decay, and fast 
cross over. We conclude in section IV, with a summary 
of our results and possible applications. We include an 
appendix, describing how to obtain analytically the scal- 
ing factor controlling the dependence of the results on 
lattice spacing. 

II. THE MODEL 

We are interested in (d+f )-dimcnsional scalar field the- 
ories with conservative dynamics defined by the Hamil- 
tonian 



(1) 



where the tree-level potential energy density is written 
for convenience as 
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We use units where h ~ c = ks ~ The parameters 
m, a, and I are positive-definite and temperature inde- 
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FIG. 1: The tree-level potential V(0) for a = 0.0, 1.5, 1.6. 

pendent. With the rescaling, </>' = (j)^/l/m, x''^ = x^m, 
and a' = a/{ra\/l)^ the potential can be written as 
V{(t>) = {m'^/l)V{4>'), with 




while the Hamiltonian in eq. [T] becomes, H[(l)\ = 
{mll)H[(t)']. Henceforth we will drop the primes. Note 
that for a — Q the model is Z2-symmetric, while for 
a — 3/2, V{4)) is a symmetric double-well. For a > 3/2 
the minimum at = becomes the false vacuum, and 
the true vacuum is at 0+ = a + \/ o? — 2. In figure [1] we 
show for several values of a. 

A. False Vacuum Decay 

Although our results could be extended to models at 
zero temperature and, thus, with false-vacuum decay 
controlled solely by quantum tunneling, we will treat 
here mostly finite-temperature models. We use a finite- 
temperature description merely for its convenience in set- 
ting the properties of the initial state, which we will take 
to be a state of thermal equilibrium prepared in the po- 
tential above with a = 0. But before we describe our 
procedure, it is convenient to briefly review some rele- 
vant results from HN. 

Consider a model with an asymmetric double-well po- 
tential, such as the one of eq. [3] with a > 3/2. The 
field is initially in thermal equilibrium with only small- 
amplitude fluctuations about the metastable minimum 
at = 0. One then computes the free energy T — 
—TlnZ, where Z is the partition function, given by 
the path integral / V[(j)\ exp[-F[(j)]/T]. At finite tem- 
perature, F[(j)\ is the Euclidean free-energy functional of 



the configuration (/>, which may contain temperature and 
quantum corrections to the tree-level theory defined in 
eq. [TJ To lowest order in perturbation theory, these cor- 
rections are included in the effective potential. At T = 0, 
one should use the {d + l)-dimensional Euclidean effec- 
tive action S[(j)], as opposed to F[(f)]/T, in the exponent 

ii- 

The calculation proceeds as follows: first, one identifies 
the stationary points of F[(j)], given by the solutions to 
the equations of motion. There are two for generic asym- 
metric double- well potentials, such as the ones in figure 
[T] the false vacuum at and the spherically-symmetric 
bounce or critical bubble (j^biT), found by solving 

dl<p+^^dA = %-, (4) 

r ocp 

where d is the number of spatial dimensions. The solu- 
tion must satisfy the boundary conditions, 06 (0) — 4>q, 
dr4>b{0) = 0, and limr^oo — 4>m- Second, adopt- 
ing the semi-classical approximation where the path in- 
tegral is dominated by small fluctuations about the clas- 
sical paths, one sums over quadratic fluctuations about 
the stationary points of the free energy using the saddle- 
point approximation to the Gaussian path integral. The 
calculation is completed by summing over the contribu- 
tion from many bubbles, which makes use of the "dilute 
gas" approximation, where the nucleating bubbles are as- 
sumed to be far enough away from each other so as not to 
overlap, being thus treated as independent events. The 
end result, the probability per volume per unit time of 
nucleating a bubble of "true" vacuum within the false 
vacuum background for a system at temperature T is, 

r ^ Tt'^+i) exp[-F[0fc]/r] . (5) 

The quotes are a reminder that the nucleating bubbles 
are not initially at the true vacuum, only close to it 
(00 7^ 0+)- The sums over bubbles and over quadratic 
fluctuations about the stationary points give rise to a 
prefactor which can be roughly approximated by T'-'*'*"^) 
or, in quantum tunneling, by Af*"*"^, where M is the rel- 
evant mass scale. Here, we will be mainly interested in 
the exponential factor which, in most scenarios, domi- 
nates the false vacuum decay rate. Summarizing, the 
calculation of false vacuum decay in HN relies on two 
key assumptions: small (quadratic) fluctuations about 
the false vacuum and about the critical bubble, and the 
dilute gas approximation. 

B. Thin Wall Approximation 

Although the equation describing the critical bubble, 
eq. m does not, in general, have exact solutions, an esti- 
mate of the energy of a critical bubble can be obtained in 
the "thin wall" approximation [l|, Q . Essentially, when 
the potential V{(t)) is nearly degenerate, the critical bub- 
ble will have a spatial extension, or "radius," {Rb) much 
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larger than the thickness of the wall separating its interior 
from its exterior, defined as the region where cj) changes 
abruptly from 0o — 0+ at r = to at r — > oo. There 
are thus two main contributions to the free energy (or Eu- 
clidean action) of the critical bubble: the wall (r ~ i?b), 
and the interior of the bubble (r < i?b), where one writes 
4>b{T) = The approximate expression for the bubble 
energy is, 



CdR 



(6) 



where |A1/| = \V{(t>+)-V{<t)ra)\: and Cd = 2ti'^/^ /T{d/2) 
is related to the volume of a d-dimensional sphere of ra- 
dius RhyVd = Cd^. Also, a = J dr [i(</.;,)2 + is 
the surface tension. Minimizing with respect to R gives: 
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and the nucleation energy-barrier in the thin wall approx- 
imation is. 
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(8) 



Within the thin-wall approximation, the bounce so- 
lution may be approximately parameterized by 4){r) ~ 
^(1 — tanh(i(r — Rb)))- Using a scahng argument, we 
can rewrite the surface tension as ct = J dr{dr(j>b)^ , which 
allows us to integrate a to obtain 



lim 

-Rb>i 
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The volume contribution from the potential is only a 
function of the asymmetry Ay = i — — a-^ — 2 + 
^(a^ -|- a'^vc? — 2). It is convenient to write a ~ 
3/2-1- da, where ^ 1 in the thin wall approxima- 
tion. Then, expanding the expressions for the energy, eq. 
[SI and the radius, eq. [T] about 5a = 0, we can obtain 
approximate analytical expressions for both as follows: 
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(10) 
(11) 



Although these expressions breakdown quite quickly as 
6a increases, they provide the dominant scaling of energy 
and radius with the coupling a. 



C. Oscillons: A Brief Review 

Alongside with the critical bubble or bounce, oscillons 
will play a key role in the mechanism for fast vacuum 
decay. In fact, the title of this work alludes to their 
presence in the metastable minimum and their effect on 
the decay rate. As such, it is useful to briefly review 
their properties. In their simplest form, oscillons are 
spatially-extended, very long-lived, time-dependent so- 
lutions of the nonlinear Klein-Gordon equation "o*], or 
of other PDEs with amplitude-dependent nonhnearities 
[T^ . Their properties have been extensively studied dur- 
ing the past decade in two [ll| , three , and higher [l^ 
spatial dimensions, and, more recently, in U(l) Abelian- 
Higgs models [3 and in SU(2)xU(l) models ^l- How- 
ever, it is fair to say that a more fundamental under- 
standing of their existence and longevity is still lacking. 

In the context of relativistic scalar field models, os- 
cillons are characterized by large-amplitude oscillations 
about the vacuum state. Assuming spherical symmetry, 
it has been shown T3] that they only exist if these fluc- 
tuations probe beyond the inflection point [^inf] of V{4>), 
although they may be more general . The reader may 
consult the references cited for more details. 

Scalar field oscillons have been found by two meth- 
ods. The first, and simpler, method makes use of an 
initial condition 0(0, r) that resembles the oscillon solu- 
tion, for example, 4>{0,r) = (/>o exp[— r^/i?^]^ fQj. y((^) of 
eq. [21 For values of i? > i?inin, and (j>o > 0inf, the field 
will evolve into the oscillon configuration, where it will 
remain for a lifetime that is sensitive to the values of 
00, ^min, and d. i?min Can be estimated analytically for 
degenerate and non-degenerate polynomial potentials for 
an arbitrary number of spatial dimensions [13| . 

The second method, more general, shows that oscil- 
lons emerge, under very general conditions, when large- 
amplitude fluctuations about a given vacuum state with 
sufficient thermal (or quantum noise) are present. In the 
context of 2d models, Gleiser and Howell showed that 
they emerge after a quench from a single- well to a double- 
well potential [13] ■ In the model of eq. [31 this would 
correspond to a change from a = 0toa=:1.5. 

Inspecting figure [1] one can see qualitatively that the 
quench induces coherent oscillations of the field's zero 
mode, 4>{t) = (1/^) / d'^x(j){x,t) about the metastable 
minimum at </) = due to a "widening" of the potential 
there. Other mechanisms may induce large-amplitude 
oscillations of the zero-mode, producing similar results. 
This happens, for example, in the context of reheating 
in certain models of cosmological inflation [l^. An- 
other possibility is a symmetry-breaking interaction that 
changes the shape of the potential from single to double- 
well, as in the case of the quench above, or one that 
simply "widens" the curvature of the potential about the 
vacuum. Quantitatively, the energy from these coherent 
oscillations is transfered to higher k modes via paramet- 
ric amplification and trigger the emergence of oscillons. 
We refer the reader to reference 17[ for details. Not 
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FIG. 2: The volume-averaged field, or order parameter, 
{(/>(t)) = <j){t) which corresponds to figO 



surprisingly, the same behavior ensues in 3c?: a quench 
induces oscillations of the field's zero mode, as can be 
seen in figure [D The results were obtained in a cubic 
lattice of volume = vifith periodic boundary con- 
ditions, lattice spacing dx — 0.5, and a changed from 
a = to a = 1.5. The initial thermal state was prepared 
at T = 0.29. 

In figure [31 we show snapshots of the field, display- 
ing the appearance of bubble-like oscillons. Shown are 
the isosurfaces at <j) — {0.75(blue), l.O(red), 1.5(green)} 
[colors online only]. The long-lived, localized oscillations 
which emerge in synchrony are the oscillon configura- 
tions. 



D. The Hartree Potential and Lattice 
Implementation 

As mentioned above, the initial state is prepared as 
a thermal state centered at in the potential of 
eq. [3] with a — 0. The thermal fluctuations will induce 
changes in the potential for (j), which, in leading order 
in perturbation theory, may be approximated by the Ho- 
mogeneous Hartree Approximation (HHA). An excellent 
account of the HHA, using successive truncations within 
the hierarchy of correlation functions to obtain corrected 
equations of motion, can be found in ref. [2(| . 

The HHA assumes that the fluctuations of the field and 
its associated momentum remain Gaussian throughout 
its evolution. Thus, the approximation works well just 
prior and just after the quench for all temperatures, and 
for low temperatures at all times. The Hartree potential 
can be derived as follows: write the field a,s (j) = (j) + 5(j). 
Then, Vh{4>) = {V{4' + S4>))j where () means an average 
over all fluctuations S(j>, with {6(j)) = and (Scj)'^) = (3. 



{5(f>^) is the translationally-invariant mean-square vari- 
ance of the field, which, in thermal equilibrium (the ini- 
tial state), is constant and proportional to T. We then 
obtain (suppressing the bar), 



(12) 



yjf((/.) = -a/3</.+ 1 + -/3 



In thermal equilibrium, 
modes in fc-space satisfy 



the momentum and field 



{\m\') 



T 



fc2 



(13) 
(14) 



where the Hartree mass is defined as = ^if(O). 

Due to the dependence of {\(f){k)\^) on the wave num- 
ber k, a lattice implementation of the field theory will 
introduce a lattice spacing {5x) dependence [2l[. We 
thus write, for an UV cutoff A, /3 = a^/T, where the 
constant a^di which changes in different spatial dimen- 
sions, can be obtained numerically by comparing the lat- 
tice results with eq. [141 It can also be obtained ana- 
lytically, as shown in the Appendix. For example, for 
6x — 0.2, we obtain, through this method, a2d = 0.518 
and a^d = 1.194, in two and three spatial dimensions, 
respectively. These values agree very well with the nu- 
merical estimate. (See Appendix for details.) 

In figure [H we plot on the left the ensemble- 
averaged two-point correlation functions for the field 
(green squares, bottom) and its related momentum (blue 
squares, top). The black lines are the analytical formu- 
las of eq. [131 On the right, we take the logarithm of 
the initial data set for the two-point correlation function 
for the field (green squares and black line), and com- 
pare to its value during the simulation. The red squares 
are the data at every half-second for the time interval 
360 < t < 400, when (f>{t) is well settled in the potential 
minimum (check figure [2]). The blue squares correspond 
to the same data, but after applying a thermal filter. 
One can clearly distinguish two populations of modes: 
low-fc, oscillon-related modes far from equilibrium (for 
k < 1.5), and modes that remain in thermal equilibrium 
for all times (fc > 1.5). 

Introducing (3 allows us to relate models with different 
temperatures and lattice spacings with a single parame- 
ter. From now on, we will refer to (3 as the effective lattice 
temperature. In figure[5l we plot Vh{4') for various values 
of (3 and a = 1.6. The interpretation of /3 as an effective 
temperature should be apparent. 

The initial state was prepared by thermalizing the field 
using standard Langevin techniques [2^ 



-I- 70 - V^(f> = -dV/d(t) + C 



(15) 



where C is a Markovian noise with two-point correlation 
function obeying the fluctuation-dissipation relation. 



(C(x,t)C(x',t')) = 27r5(x-x')<5(t-t') 



(16) 
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FIG. 3: Snapshots of a 3ci scalar field simulation after a quench from a single to a symmetric double well potential (a = — > 1.5). 
The oscillatory bubble-like configurations which emerge in synchrony are identified as oscillons. The lattice parameters were 
dx = 0.5, T — 0.29 and L"^ — 64"^. Plotted is the = {0.75, 1.0, 1.5} isosurfaces in blue, red and green, respectively. The 
plots are snapshots at t = {45.5,51.5,54.5,60.5,62.0,194.5}, respectively, increasing left to right and top to bottom. Note 
that oscillons can be tracked over a period of oscillation as they reappear in approximately the same place about ~ 9s later. 
For example, the top right corner of the first, second, fourth and fifth slices show the same oscillon. The last snapshot shows 
oscillons still present at late times. They eventually disappear as thermal equilibrium is restored. 



The initial state is solely determined by the tempera- 
ture of the bath T. Note that this thermalization is sim- 
ply a device to set an initial distribution for the field 
modes satisfying eq. [M] . One could have used alterna- 
tive T-independent methods to set up initial conditions. 
Ours is in tune with the more conventional statement 
that vacuum decay initiates from a metastable thermal 
state. This initial thermalization was verified by com- 
pliance with equipartition, such that the average kinetic 
energy of each mode was T/2. 



We used a standard leapfrog algorithm [23| on cubic 
lattices with volume V — — (256 * dx)^ with periodic 
boundary conditions. We verified that our results do not 
show any relevant dependence on L. The lattice was 
evolved in time-steps of 5t = 0.025. We checked that the 
Hamiltonian evolution conserves energy to 0{5t^). 



III. POWER-LAW VACUUM DECAY 

Within the framework of HN, the nucleation energy 
barrier can be calculated numerically by integrating the 
equation: 

a20 + fcila,0 = a^^ff(^,T). (17) 

In figure[6l we plot the energy barrier Eb{a, [3) as a func- 
tion of asymmetry a for several values of (3 for d — 2>. For 
reference, we also included the Euclidean bounce action 
for quantum decay at T = 0. It is apparent that the 
different curves are related by a /3-dependent scaling, 

Eb{a,l3)^Eb{a,Q)A{j3), (18) 

where we found that A{P) = (1 — cq;/3)^, with c = 2.75 
giving the best fit. In order to test this approximation, in 
figure [3 we plot the ratio Eb{a,P)/ [Eb{a,0)A{(3)]. It is 
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FIG. 4: The radially-averaged 2-point field and momentum correlation functions. On the left, we show the results for the 
thermal initial conditions. The black continuous lines correspond to (|7f(fc)p) = T{Sx)~^ (top) and to (|(^(fc)p) — T5x~^ /{k^ + 
mn} + co\/k (bottom). (The small constant co = 0.024 used to fit the analytic prediction with the lattice results is due to the 
abrupt lattice cutoff introduced by 5x.) The blue and green squares are from the simulation. On the right, we plot the logarithm 
of the two-point correlation function for the field. The black line and green squares correspond to the thermal initial state plot 
on the left. The red squares are the data plotted at every half second for the time interval 360 < t < 400. The blue squares 
are also the data, after applying a thermal filter. One clearly distinguishes two populations of modes: the oscillon-related, low 
k modes (fc < 1.5), which sharply depart from the thermal spectrum, and the purely thermal modes, which remain in thermal 
equilibrium throughout the simulation {k > 1.5). These figures were generated using the data from the same run as in figlS] 
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FIG. 5: Hartree potential for a = 1.6 and several values of 
/3 = 0,0.1,0.2,0.4,0.75. 



clear that the scaling holds to within 10% for the range 
of parameters where resonant nucleation ensues. The 
scaling allows us to relate the finite-temperature bounce 
action computed with the Hartree approximation of eq. 
1121 with the tree-level bounce action computed with eq. 



FIG. 6: The logarithmic energy of a 3d bounce or critical 
bubble, ln[Eb(a, 13)], for (3 = {0.0,0.025,0.05,0.10,0.125} as 
a function of a. There is a clear scaling with /3 that worsens 
as (3 is increased. For /3 = 0, the bounce "energy" is the d — 4 
Euclidean action. 



m 



7 




CO 1.5 - 




FIG. 7: Testing the scaling of the bounce energy with /3: we 
plot the ratio Eb{a, P)/ [Ei,{a,Q)A{/3)] versus a for various 
values of /3. The scaling holds well (dashed lines denote 10% 
accuracy) and worsens as l3 increases. 



Computing the Exponent in Resonant 
Nucleation 



According to HN theory, the decay rate is given by eq. 
[SI Is this also true after an abrupt change in the potential 
V{(f>) - promoted either by a quench, a new interaction, 
or by a coupling to another field that induces oscillations 
on (f>{t), as in certain inflationary models [H, [l^? In ref. 
[2^ (GH), it was found that, in 2d, changes that induce 
large-amplitude fluctuations of the field's zero mode may 
drastically alter the decay rate from an exponential to 
a power law for a range of asymmetries (a) and tem- 
peratures (T). In GH, it was argued that the decay is 
triggered by the resonant nucleation of oscillons, which 
may coalesce to form a critical bubble or that simply 
become unstable and grow to become a critical bubble 
themselves. This mechanism was dubbed resonant nucle- 
ation (RN) , and was found to lead to a power-law decay 
per unit area 



FIG. 8: The value of the decay power _B as a function of /? in 
two dimensions. 



r - (Eb/Ty 



(19) 



with B ~ 2.464±0.134 for the range of parameters where 
RN is operative. In figure[5]we display our results for the 
exponent B m d = 2. The lattice spacing was Sx = 0.2, 
and the lattice size 2048^. Results were averaged over 
10 runs. Errors were smaller than boxes. [Note that this 
value of B and the error bars are smaller than those 
quoted in Also, we found that B is only weakly de- 
pendent on p. This discrepancy is due to better statistics 
and to a more precise method of distinguishing between 
RN and cross-over transitions, to be explained below.] 

In order to test if a power-law behavior holds in 3d, 
we repeated the GH procedure for different values of the 
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FIG. 9; The volume- averaged field, or order parameter, (j}{t) 
for Q = 1.589 1.78 for /3 = 0.085. [Results are better seen 
in color. 1 



asymmetry a and the lattice temperature parameter /3. 
The field is initially prepared in a thermal state in a 
single well (a = 0). The potential is then changed to 
an asymmetric double well, by setting a > 3/2. The 
subsequent dynamics of the field is conservative, that is, 
we set 7 = in eq. 1151 We then measured the volume- 
averaged value of the field, 4>{t)- Results for (3 = 0.085 
and several values of the asymmetry a are shown in figure 

M 

One can easily see that, for values of a close to degener- 
acy (towards the right of the figure), (t>{t) performs many 
damped oscillations about the metastable minimum of 
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Vni't') before irreversibly transitioning to the gfobal min- 
imum. For clarity, we cut olf the plot before the decay 
completes. As the asymmetry is increased (towards the 
left of the figure) , the number of oscillations decreases un- 
til the transition completes by cross-over, that is, within 
a time-scale which is shorter than a typical oscillation 
period. [More on this soon.] 

Adopting as the decay time-scale (r) the time at which 
(j){t) crosses the maximum of VH{<t>), we can plot the loga- 
rithm of decay times as a function of the logarithm of the 
bounce energy {Ei,/(3), where Eb is computed from the 
solution of eq. [17] for a given pair (a, /3). The results are 
ensemble-averaged over 50 runs, and displayed as indi- 
vidual points (diamonds) in the left plot of figure [TUl On 
the right side, we plot the logarithm of {Eb/P^) to show 
the simple scaling with /?. For reference, if Sx = 0.25, 
the value we adopted in our simulations, the tempera- 
ture would be in the range 0.09 < T < 0.14. [To get 
the value of the temperature, use that T = (S/a^d- An- 
alytically, a^d = 0.955 (see Appendix), and numerically 
a-zd — 0.933, an error of only 2.4%.] The slopes of the 
straight-line fits are the numerical values for the decay 
power B. 

Note that we are not including decay times smaller 
than r ~ 22. Times shorter than this are characteristic 
of a cross-over transition and would not fall under reso- 
nant nucleation. This can be made clear by investigating 
how we extracted the slopes of the log-log plot, that is, 
the constant B controlling resonant nucleation. As an 
illustration, fixing (5 = 0.0954, consider the results for 
the decay times displayed in figure [11] From the data, it 
is clear that there is a sharp departure from a straight- 
line fit below ln(T) ~ 3.2. This happens consistently for 
all parameters we investigated, offering a natural cutoff 
below which the transition occurs via cross-over. 

In figure [T^ we show our results for the power B con- 
trolling resonant decay. We have included two data sets: 
the squares are the results obtained using the Hartree 
effective potential to compute the bounce action and to 
read off the nucleation time-scale (which depends on the 
location of the maximum). The circles are the results 
obtained using the tree-level potential. It is clear that 
for Vff ((/)), B is fairly independent of (3 within the range 
where resonant nucleation applies. The average value is 
B = 1.327 ± 0.059. The values agree well for (3 > 0.09. 
For smaller /3, it is necessary to increase the value of a 
to probe into the resonance nucleation regime, leading to 
larger errors in the tree-level estimates. 

We also verified that the value of B is not very sensitive 
to changes in lattice spacing. As an illustration, let us 
choose Sx = 1. For T ~ 0.55, which gives f3 — 0.131, we 
repeated the numerical experiments to obtain B ~ 1.68. 
By adjusting T, this value can, for example, be compared 
with the one obtained with lattice spacing Sx = 0.25 for 
the same (3 = 0.131, which we measured to be B = 1.47. 
This is a change of only 13% in B for a factor 4 change in 
lattice spacing. As a second illustration, we compared the 
results for Sx = 0.35 and T = 0.175, for which /3 = 0.119, 



with Sx = 0.25, at same (3. The values for the power-law 
decay were B — 1.458 and B = 1.353 respectively, a 7% 
difference. Within the accuracy of our simulations, we 
conclude that B is very weakly dependent on Sx and, 
most importantly, on f3. 

B. Phase Diagram for Vacuum Decay 

We summarize our results in a phase diagram display- 
ing the three possible mechanics for false vacuum decay 
as a function of the parameters controlling the asymme- 
try of the potential (a) and the temperature of the ini- 
tial metastable state (/3). Homogeneous nucleation (dark 
blue region labeled as HN) describes vacuum decay for 
low enough temperatures for all values of asymmetries. 

Resonant nucleation (RN) falls within the region 
bracketed by the two left to right curves (light blue re- 
gion labeled as RN). The diamonds and circles are the 
results of the simulations. The circles, being close to the 
cross-over region (r 25), were not used to compute the 
decay power B. The curves bracketing RN from above 
(cross over, purple region) and below (HN region) were 
obtained by using the average value of the field, 0(<), as 
a guideline: a field well-localized within the false vacuum 
{(f>{t) <C (pini) will only decay via HN; a field for which 
(f>{t) ^ (pini will easily cross over. RN lies in between. 
Noting that the amplitude of thermal fluctuations about 
(i>{t) is given by \\/JS(j^\ = [•\//3/2|, we obtain the con- 
dition, 

VHi^inf±y/m) = 0. (20) 

This condition translates into a relation between a and /3 
that generates the two curves bracketing the RN region 
horizontally. 

The steeper line to the left is obtained by examining 
the dynamics of resonant nucleation. As noted in ref. 
p^ ]. in the region near degeneracy, where the critical 
bounce is large, RN will occur due to the coalescence of 
two or more oscillons. In figure 1141 we compare the ra- 
dius of a bounce with that of oscillons as a function of a 
for f3 = 0.107. There are three distinct regions: For large 
a ^ 1.7, the bounce and oscillon practically match. This 
region is typical of quick cross-overs or at most border- 
line RN, where a single oscillon becomes unstable, grows, 
and becomes critical. In the intermediate region, where 
1.56 < a < 1.7, i?b ounce 

< 2i?osc: here, two oscillons may 
coalesce to form a critical bubble. This is the region typi- 
cal of RN. Lastly, for a < 1.56, i?bouncc > 2i?osc and more 
than two oscillons are needed to produce a critical bub- 
ble, a process that, although not impossible, is strongly 
suppressed in the allowed time scales where RN can act, 
that is, while oscillons are present in the system, typically 
while dt<f>{t) 7^ 0. In this region, HN is the favored decay 
route. The steeper line in the phase diagram is thus ob- 
tained by tracing the points where i?bouncc — 2i?osc for 
each /?, in excellent agreement with the simulations: RN, 
denoted by the diamonds, ends on this line. 
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FIG. 10: Decay time-scales in 3d: we plot ln(r) vs. ln{Eb/f3) [left], and ln{Eb/f3'') [right]. From top to bottom on the left 
figure, /? = {0.084, 0.095, 0.101, 0.107, 0.119, 0.131}. 
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0.14 



FIG. 11: Extracting the decay exponent B: shown are the 
raw data points for the decay time Inr vs. (InEt/P) with 
l3 — 0.0954. We ran 50 experiments for each value of a (or, 
equivalently, of Eb). In this illustration, we used the tree-level 
potential of eq. |3]to compute Et- The straight-line slope fit 
gives B = 1.473. Note the clear departure from a straight- 
line fit for Inr ^ 3.2. We characterize time-scales faster than 
those as typical of cross-over transitions. 



In figure [TH we show the coalescence of oscillons lead- 
ing to a critical bubble that grows to complete the vac- 
uum decay. This figure should be contrasted with fig- 
ure [3] for a symmetric double- well potential. The regions 
displayed are isosurfaces at (/> = {0.9,1.75,1.85} which 
are colored blue, red, and green, respectively. We used 
a = 1.545 and (3 = 0.134. In figure [TBI we show the cor- 



FIG. 12: Results for the decay power B. Slopes as measured 
individually, as in figure 1111 The circles are calculated with 
the tree-level potential and the squares with the Hartree effec- 
tive potential. Error bars from ensemble-averaging fit within 
the sizes of the squares and circles. 

responding evolution of the volume-averaged field, 
The arrows denote the locations of the isocurvatures de- 
picted in figure [TSl 



IV. SUMMARY AND OUTLOOK 

We have presented a detailed numerical study of the 
decay of metastable vacua in scalar field theories. In 
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FIG. 13: Phase diagram for false vacuum decay. Shown are 
the three regions characterizing vacuum decay as a function 
of the asymmetry in the Hartree potential a and the effective 
lattice temperature /3. HN (dark blue) represents decay by 
homogeneous nucleation of critical bubbles. RN (light blue) 
represents decay by resonant nucleation, either by the coales- 
cence of two oscillons (diamonds) or by the unstable growth 
of a single oscillon (circles). The top region represents fast 
cross-over, within time-scales of r < 25. [The regions were 
labeled for clarity with B&W display. 




1.5 1.6 1.7 1.8 1.9 2 

CI 

FIG. 14: Comparing the sizes of oscillons [J?osc - continuous 
(blue) lines] and critical bubbles [/Jbounce - dashed (red) line] 
as a function of asymmetry a: 2iiosc is the upper line, and 
i?osc lower continuous (blue) line. Data obtained for /3 = 0.107 



particular, we investigated the effects of large-amplitude 
fluctuations about the vacuum state on the decay rate. 
These large-amplitude fluctuations were produced by a 
rapid quench from a single to an asymmetric double well 
potential. We have argued that they are characterized by 
spatially-extended, long-lived oscillatory configurations 
called oscillons. There are, of course, short-lived, small- 
amplitude configurations as well, but those play a small 
role on the dynamics of vacuum decay. They are incor- 
porated, via a perturbative approach, into the effective 
Hartree potential. 

Our main result was to obtain a phase diagram charac- 
terizing the three possible mechanisms for vacuum decay. 
Two of them, slow homogeneous nucleation and rapid 
cross-over, are well known. Homogeneous nucleation en- 
sues when fluctuations about the metastable vacuum are 
small-amplitude, within the perturbative regime. Cross- 
over ensues when fluctuations are very large-amplitude, 
probing well beyond the inflection point of the effec- 
tive potential. The novel mechanism, resonant nucle- 
ation, lies between these two and ensues when fluctua- 
tions probe the neighborhood of the inflection point. We 
have shown that, as in 2d [l^], 3d resonant nucleation is 
characterized by a power-law decay rate, r ~ {Eb/(3)^ , 
where i?b is the energy of the critical bubble computed 
with the effective potential, (3 is the effective lattice tem- 
perature, and B = 1.327 ± 0.059 is the exponent con- 
trolling the decay. We argued that resonant nucleation 
typically occurs when two oscillons coalesce to form a 
critical bubble or, marginally, when a single oscillon be- 
comes unstable to growth. Although the coalescence of 
three or more oscillons is possible, this process tends to 
be highly suppressed due to its long time-scale. 

Our results open many avenues for future investiga- 
tion. Can oscillons play a role in inflationary cosmol- 
ogy? During reheating, the inflaton undergoes damped, 
large-amplitude oscillations about the potential mini- 
mum. This is essentially the same mechanism that gives 
rise to oscillons. Of course, one must add the effects of 
the cosmic expansion, but preliminary results in simple 
Id models indicate that oscillons not only will be present 
but will actually be stabilized by the expansion, becom- 
ing spikes in the energy-density distribution [2^ . We are 
currently investigating such possibility in a full 3d simu- 
lation. Another area where the effects of oscillons should 
be further investigated is during spontaneous symmetry 
breaking. In the present work, we investigated a sim- 
ple Z2 symmetry breaking, as we quenched from a single 
to an asymmetric double well potential. It is clear that 
the dynamics of symmetry breaking is very sensitive to 
the model parameters and that oscillons can play a key 
role. In a related work, we found that low-momentum 
vortex-antivortex scattering in broken U(l) gauge mod- 
els can lead to very long-lived oscillons characterized by 
a (practically) non-radiating oscillatin g rn agnetic dipole 
and nontrivial Chern-Simons number [1^. Farhi et al. 
and Graham 15] found a long-lived oscillon in broken 
SU(2)xU(l) models when the Higgs mass is twice the 
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FIG. 15: Snapshots of a 3d scalar field simulation after a quench from a single to an asymmetric double well potential {a — 1.545, 
T — 0.28, dx = 0.5, l3 = 0.134). The oscillons coalesce to become a critical bubble that grows to complete the vacuum decay. 
The slices are at times t — {27.0, 38.0, 45.5, 54.5, 71.0, 79.5}. The blue, red and green isosurfaces are at = {0.9, 1.75, 1.85}, 
respectively. The average value of the field along with the respective isosurface locations are also noted on figfT6l 



mass. If this turns out to be the right mass ratio, as hope- 
fully we wiU soon know from the LHC, we should expect 
electroweak oscillons to exist in Nature. Alternatively, 
they may exist even for a wider range of values. Taken 
together, these results demonstrate the rich physics of 
time-dependent, spatially-extended field configurations. 
This richness is only beginning to be explored. 
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V. APPENDIX: ANALYTICAL ESTIMATE OF 
EFFECTIVE LATTICE TEMPERATURE 

Consider the simple case of a quadratic potential 
Vq — , for which the Hartree potential is simple 
Vh — + \P- Consider also the 1-loop potential with 
UV cutoff A, 

^il(0) = Vo + ^£ ^In {p' + V^') . (21) 
In d = 2 the integral gives, 

^il(0) = Vo + {l - In (^) } , (22) 

while in d = 3 wc obtain 

^il(0) ^Vo + ^ {w^'A - 7r(yo")'/'} ■ (23) 

Focusing on the simple quadratic potential, for which 
Ig" = 1, we now match the cutoff-dependent term of Vil 
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lime 

FIG. 16: The volume- averaged field, or order parameter, 
{(j>{t)) = cj){t) for a = 1.545 for /3 = 0.1385. The colored 
arrows label the locations of the isosurfaces depicted in figure 

(multiplied by d/2) with the /^-dependent term of Vh (eq. 
to obtain, 

/3=|(l + ln(|^)).a..T, [2d] (24) 

and 

ST 

f^=7ZIZ^ «3.r, [3d] (25) 



where we used A = Tr/{Sx) and introduced, for conve- 
nience, the lattice-dependent constant a^d- For example, 
for Sx — 0.2, a2d = 0.518 and 03^ = 1.194. This value 
agrees very well with the numerical estimate, obtained 
by comparison with eq. [TH 

After comparing terms, one can write most gener- 
ally as 



where ~ implies take the first term of the Taylor series 
expansion of the integral should be a^d- In the Table 
below, we compare the numerical and analytical (eqn. 
I26|) measurements of a^d in different dimensions. 





analytical 


numerical 


lattice size 


aid 


0.5 


0.5 ±0.001 


10^ 


a2d 


0.518 


0.53 ±0.005 


10242 


asd 


1.194 


1.19 ±0.005 


100^ 


O-id 


3.125 


3.84 ±0.005 


32* 



With this simple analytical method, it is possible to 
relate the initial state temperature T to the effective 
lattice temperature (3 for any choice of lattice spacing. 
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